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THE  UNSTEADY  FLOW  WITHIN  A  SPINNING  CYLINDER 

1.  ABSTRACT 

A  theoretical  analysis  is  given  of  the  unsteady  flow  within  a  liquid-filled 
cylinder  of  finite  length,  which  is  started  impulsively  to  spin  about  its  axis. 

It  is  established  that  a  secondary  flow,  which  is  caused  by  the  cylinder 
ends,  has  a  remarkable  effect  on  the  generation  of  spin  in  the  liquid. 

The  theoretical  results  are  compared  with  experimental  observations  and 
good  agreement  is  found, 
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I.  INTRODUCTION 


It  is  veil  known  that  spin- stabilized  shells  can  become  dynamo cally  un¬ 
stable  if  they  are  filled  with  a  liquid.  According  to  a  theoretical  analysis 
of  K.  Stewart son  (1959)>  the  stability  of  a  shell  containing  a  cylindrical 
liquid-filled  cavity  can  be  predicted  if  the  liquid  is  in  rigid  body  rotation. 


However,  for  a  liquid  of  small  viscosity,  it  requires  a  relatively  long 
time  for  the  liquid  to  reach  full  spin,  and,  during  the  transition  period,  the 
shell  might  become  dynamically  unstable  even  though  it  might  be  stable  at  its 
final  state  if  the  liquid  attained  rigid  body  rotation.  In  the  course  of 
experimental  investigations  (E.  G.  Karpov,  1962),  severe  dynamic  instabilities 
of  liquid-filled  spinning  shells  have  been  observed  i-n  cases  where  the  shell 
should  have  been  stable  according  to  Stewartson's  theory  and  the  assumption  of 
rigid  rotation  of  the  liquid  filler.  In  order  to  extend  the  prediction  of  in¬ 
stabilities  in  such  cases,  it  appeared  desirable  to  analyze  the  problem  of  un¬ 
steady  fluid  motion  within  a  cylinder,  started  to  spin  about  its  axis. 


If  the  cylinder  is  infinitely  long,  a  solution  to  the  problem  is  obtained 
without  difficulty,  but  the  expectation  that  this  solution  might  approximately 
be  valid  for  slender  but  finite  cylinders  proves  to  be  wrong.  It  is  found  that 
the  effect  of  the  cylinder  ends  on  the  fluid  motion  is  not  only  not  negligible, 
hue  dominating.  The  fluid  motion  is  entirely  changed  by  the  presence  of  a 
secondary  flow  Induced  by  the  cylinder  ends.  The  secondary  r'inw  convects 
spinning  fluid  from  the  walls  into  the  inter ‘or  of  the  cylinder  and,  as  a 
consequence,  the  fluid  attains  rotational  motion  many  times  faster  than  without  % 
secondary  flow.  In  the  following,  a  theoretical  analysis  is  given  of  the  un- 

J 

ctcadj  flow  within  a  cylinder  which  is  started  +o  spin  about  its  axis  of  rotation! 
The  results  are  then  compared  with  some  experimental  data  obtained:  l)  from 
spin  decay  data  of  liquid-filled  shells;  and  2)  from  direct  observations  of  the 
secondary  flow  within  a  spinning  transparent  cylinder. 


II.  THEORETICAL  ANALYSIS  OF  THE  FLOW 


a .  Structure  of  the  Secondary  Flow 

The  diagram  in  Figure  1  shows  an  axial  section  and  a  cross-section  of  the 
spinning  cylinder.  The  height  of  the  cylinder  is  h,  the  radius  a.  For  the 
.'or rowing  analysis,  we  use  a  polar-coordinate  system  6,  r,  z  with  the  origin  in 
the  center  as  shown  in  the  diagram  (Figure  l).  The  velocity  components  are  v,  u, 


respectively.  The  Navier-Stokes  equations  in  these  coordinates  for  a  flow  of 


rotational  symmetry  are: 


dv 

/dv 
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(lb) 

(lc) 

(ld) 


The  boundary  conditions  are  u  =  w  =  0,.  v  =  ro>  at  z  =  -  h/2,  and  u  =  w  =  0  v- 

at  r  =  a,  with  u  =  0  at  r  =  0, 


Let  us  assume  the  cylinder  is  started  at  time  t  =  0  .to  spin  about  its  axis 
of  rotation  with  the  constant  or  time-dependent  angular  velocity  <n. 

If  the  cylinder  is  infinitely  long;  i.e.,  h  ->  °°,  the  Equations  (lb),  (lc)  and 
(Id)  with  the  eorresponding  boundary  conditions  are  solved  by 


and  Equation  (la) reduces  to  a  linear  differential  equation: 

2 

dv  3  v  ,  dv/r  /nN 

3t  =  v  £S  +  -ST  '  (2) 

where  v  depends  on  r  and  t  only;  i.e.,  v(r,t).  No  such  solution  with  u  =' w  =  0 
is  possible  if  the  cylinder  has  a  finite  length .  In  the  vicinity  of  the  endwallB, 
z  ■  +  h/2,  the  circumferential  component  v  must  depend  on  z  and,  according  to 
Equations  (lb)  and(lc),  the  velocity  components  u,  w  must  be  different  from  zero. 

In  fact,  the  fluid  particles  at  the  cylinder  ends  rotate  with  the  velocity  of 
the  walls  and,  therefore,  are  subject  to  centrifugal  forces.  Due  to  these 
centrifugal  forces,  the  particles  close  to  the  endwalls  are  driven  outward, 
creating  a  secondary  flow  with  velocity  components  u,  w.  We  can  assume,  however, 
with  the  reservation  of  a  final  proof,  that  the  secondary  motion  u,  w  is  very 
slow  except  ba  a  thin  boxmdary  layer  region  at  each  of  the  endwalls .  Thus,  we 
can  divide  the  entire  flow  region  in  two  parts,  the  boundary  layer  region  close 
to  the  endwalls  and  the  rest,  of  the  flow,  which  we  will  call  the  core  flow.  A 
similar  flow  structure  was  found  by  Ludwieg  ( 19!? l)  for  the  steady  flow  in  a 
rotating  duct,  and  the  following  considerations  will  be  analogous  to  those  given 
by  Ludwieg  (l95l). 


■am 


a 


b.  Approximate  Equation  for  the  Core  Flow 

Let  us  use  the  notation  vq,  u  ,  wQ,  pQ  for  the  velocity  components  and  the 
pressure  in  the  core  flow.  If  we  apply  Equation  (lb)  to  the  core  flow,  we  can 
neglect,  the  terms  containing  u,  w,  according  to  our  assumption  that  the 
secondary  motion  is  very  slow  in  the  core  flow,  and  we  have  approximately 


(3) 


From  Equation  (lc)  we  observe  that  the  axial  pressure  gradient  is  very 
small;  or,  within  our  approximation,  we  can  assume  that  the  pressure  is  independ 
ent  of  z.  But  then,  according  to  Equation  (5),  vq  must  be  independent  of  z,  and 
Equation  (la)  applied  to  the  core  flow  gives: 


(*0 


Since  v  is  independent  of  2,  it  can  be  seen  from  Equation  (4)  that  also  uq  must 
be  independent  of  z.  Thus,  Equation  (4)  reduces  to  a  partial  differential 
equation  in  the  independent  variables  r  and  t. 

Before  we  can  solve  Equation  ( 4)  for  the  circumferential  component  vQ,  we 

must  have  an  additional  delation  which  allows  us  to  express  the  other  dependent 

variable  u  by  v  .  The  required  additional  relation  between  u  and  v  will  be 
0  0  0  0 

given  by  the  coupling  between  boundary  layer  flow  and  core  flow. 


o.  Boundary  Layer  Flow 

Since  the  flow  within  the  spinning  cylinder  must*  be  symmetric  with  respect 
to  the  middle  plane  z  =  0,  we_ean  restrict,  our  analysis  to  the  boundary  layer  at 
one  of  the  endwalls,  say  z —  h/2. 


The  boundary  layer  equations  can  be  obtained  from  the  Navier-Stohes 
Equations  (la)  to  (id),  applying  the  usual  boundary  layer  simplifications.  The 
radial  pressure  gradient  within  the  boundary  layer  can  be  replaced  by  the 


pressure  gradient  of  the  "outer"  flow;  i.e,,  the  core  flow  in  our  case 
according  to  Equation  (5),  ^ 


=  v?/r. 
o' 


Thus, 

^v 

The  friction  forces  reduce  to  v  — 5  > 


etc.  Although  the  boundary  layer  flow  is  unsteady,  we  can  treat  it  as  a  quasi- 
steady  flow;  i,e.,  the  local  acceleration  terms  ,  etc.,  can  be  neglected. 
Apart  from  a  very  short  acceleration  period,  after  the  cylinder  is  started  to 
spin,  the  local  acceleration  termB  are  very  email  c^-ipared  to  the  convective 
terms.  During  the  acceleration  period,  the  flow  at  each  of  the  endwalls  is 
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essentially  the  same  as  the  flow  on  an  impulsively  started  rotating  disk.  The 
unsteady  boundary  layer  flow  on  an  impulsively  started  rotating  disk  was  in¬ 
vestigated  by  Thiriot  (l9*t0).  According  to  Thiriot’s  solution.,  the  duration  of 

2  1 

the  acceleration  period  is  about  t  «  —  ;  i.e.,  after  a  fraction  —  of  a 

u>  it 

revolution  the  boundary  layer  flow  is  almost  steady.  We  can,  therefore,  ignore 
the  acceleration  period  and  consider  the  boundary  layer  flow  as  quasi- steady  for 
all  time.  The  boundary  layer  equations  then  are: 


du 

du 

2  v 

V  _0 

52u 

u  3F 

+  w  31  ' 

r  r 

=  V 

Sz2 

(5a) 

c*v 

.  v  .  dv 

+  u-  +  w^  =  v 

SSv 

5z2 

(5b) 

5ru 

,  drw 

>0 

ST 

+  3T  " 

(5c) 

For  convenience  we  change  our  coordinate  system,  so  that  the  lower  endwall  of 
the  cylinder  is  given  by  z  =  0.  The  boundary  conditions,  then,  are: 


v  =  rai 
u  »  0 
w  =  0 
v  =.  v 

o 

u  =  0 


(3d) 


vq  enters  into  our  boundary  layer  problem  twice;  first,  it  occurs  in  Equation 
(5a);  and  secondly,  it  enters  into  one  of  the  boundary  conditions.  For  any 
“given" out er^flowir  (r);  the  boundary  layer-flow  is  determined  by  the  Equations 
(5a,  b,  c)  and  the  boundary  conditions  (5d).  Thus,  we  have  a  coupling  between 
the  boundary  layer  flow  and  the  core  flow. 


The  boundary  layer  Equations  (5a,  b,  c)  have  been  the  subject  of  many 
investigations,  von  Karmen  (1921)  considered  the  flow  on  a  rotating  disk  in  a 
fluid  at  rest  (v  =  0)  and  obtained  approximate  solutions  using  the  integral 
method  he  invented.  A  more  accurate  solution  to  the  same  problem  was 
calculated  by  Cochran  (195*0.  'The  problem  of  rigid  body  fluid  rotation  over  a 
stationary  disk  was  solved  by  Bodewadt  (19*<0). 


Batchelor  (].95l)j  Stcwartoon  (I950)i  Rogers  and  Lance  (i960)  and  others 
Investigated  the  more  general  problem  of  a  fluid  in  rigid  body  rotation  over  a 
rotating  disk.  A  common  feature  ox  the  above  mentioned  flows  is  that  they  have 
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similarity  solutions,  where  the  velocity  components  Luke  the  form  v  =  r  G(z), 
u  =  r  F(z),  w  =  h( z ) .  These  similarity  solutions  are  also  solutions  or  the  exact 
Navier-Stokes  equations,  since  the  terms,  which  are  commonly  neglected  in  boundary 
layer  theory,  vanish  identically.  Ludwieg  (l95t)  and  Squire  (1953)  linearized 
the  boundary  layer  Equations  (5a,  b,  c)  for  the  case  of  small  disturbances  about 
a  state  of  rigid  rotation.  In  the  linearized  form,  the  Equations  (5a,  b,  c) 
reduce  to  a  set  of  ordinary  differential  equations,  which  are  linear. 


For  general  vq  -  cListritiutlons,  when  neither  linearization  nor  the 
assumption  of  similarity  flow  is  applicable,  approximate  solutions  may  be  obtain¬ 
ed  by  using  the  momentum  integral  method.  Mack  (1962)  (1963)  has  applied  the 
momentum  integral  method  (1962)  and  a  simplified  momentum  integral  method  (1963) 
to  rotating  flows  on  a  stationary  disk.  The  latter  method,  which  makes  computa¬ 
tions  easy,  could  be  extended  to  our  case  of  rotating  flows  on  a  rotating  disk. 

While  lue  momentum  Integral  method  dues  nut  give  the  exact  shape  of  the 
velocity  profiles,  it  provides  fairly  good  approximations  to  certain  integral 
values,  e.g.,  the  radial  mass  flow  within  the  boundary  layer,  which  is 


r6 

M(r)  =  2jtr  .  p  /  u(r,z)  dz  where  5  is  the  boundary  layer 

thickness.  When  the  radial  mass  flow  distribution  M(r)  has  been  determined  -  for 
a  given  distribution  of  vo(r)  -  the  radial  velocity  in  the  core  flow,  uQ(r),  can 
be  found, 


Making  use  of  the  condition,  that  the  total  radial  mass  flow  (within  the 
two  boundary  layers  and  the  core  flow)  must  be  zero,  one  obtains 


2itr  .  p|" J  u(r,z) 

o 

-  |  uo(r)  =  J  u(r, z) 


dz  -L  h  u  (r) 


(6) 


dz 


Thus,  we  have  a  functional  dependence 
uQ(r)  =  F[vQ(r)] 


(7) 


This  means,  for  any  given  distribution  of  v  (r)  we  can  find  the  distribution  of 
u0(r). 


11 


In  principle  we  could  now  express  uq  in  Equation  (1+)  by  vq  making  use  of 
the  functional  dependence  {?)-  But  aside  from  the  fact  that  we  cannot  give  an 
explicit  formula  for  uq ( r ) ,  the  relation  (7)  will  be  much  too  complex  to  enable 
us  to  solve  Equation  ( 4) .  Thus  we  have  to  confine  ourselves  to  a  simple 
approximation  of  the  relation  (7). 

d.  Approximate  Formula  for  un(r) 

At  the  beginning  of  the  fluid  motion  the  circumferential  component  of  the 
core  flow,  v  ,  is  zero  and  the  boundary  layer  problem  reduces  to  the  problem  of 
the  rotating  disk  flow,  which  was  solved  by  Cochran  (1934). 

According  to  Cochran's  solution  we  have: 

» 

J  udz  »  0.443 J •  (raj)  (8) 

-  -  O  " 

and  hence ,  from  Equation  (b) 

-  |  uQ  =  0.443^"  rcn  (9) 

'l 

If  on  the  other  hand  the  fluid  finally  attains  the  Btate  of  rigid  rotation 

(vq  =  res)  the  boundary  layer  Equations  (5)  have  the  trivial  solution  v  =  jxc, 

u  n  w  a  0  and  hence  u  =  0. 

o 

The  simplest  possible  approximation  for  general  vq  distributions  then  is 
to  esHume,  that: 

-  |  u0  “  (T“  --0-)’- _  -  .  —  C3:0) 

which  is  a  linear  interpolation  between  the  two  extreme  cases. 

Wc  can  test,  the  validity  of  this  approximation  in  a  few  other  cases.  If 

the  core  flow  is  almost  a  rigid  body  rotation  with,  the  angular  velocity  oj,  i.e., 

v  =  rai  +  v '  and  |v*  I  <  <  ru>.  tl 
o  o  I  o  I 

linearized.  In  doing  this  we  t: 

v  =  ra>  4  v '  j  u  = 
o  o 

=  rtc  +  v* ;  w  - 
wlierc  the  primed  quantities  are 
(l^n,  b)  and  neglecting  terms  of 
Equations  (la,  3b)  become: 


ie  boundary  layer  Equations  (9a,  b,  c)  can  be 
ransform: 


(11) 


small  compared  with  ran  Substituting  (ll)  into 
higher  than  first  order  in  the  primed  quantities, 
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V 


2bju» 


0 


' 

2 

dz 


■s2  t 

v  +  2tt)(v'  -  v  ')  =  0 

Jte  ° 


(-12) 


These  linearized  equations  wei'e  used  by  Ludwieg  (1951)  for  the  boundary  layer 
flow  in  a  rotating  duet. 

With  the  boundary  conditions: 


v« 

=  ol 

V 

V  at  2  =  O 

u' 

-  0 1 

u 

v'  =  v 


at  z  -  ® 


Equations  (12)  are  solved  by 


v( 


V  1  1  1  -  cos 


u'  =  -  VQ'  sln^z.  e 


Jl ».  “ ) 


(13) 


For  the  radial  component  in  the  core  flow  u  we  have,  according  to  Equation  (6) 


-I  “<,<•■>  vJ-i 


If  we  replace  v  1  according  to  (ll)  by  -  (ru)  -  v  ),  we  have: 


-  7*  u  =  O.'jOO  .  (rto  -  v  )  (1*0 

2  o  no>  o 

'Hi is  formula  for  uq  is  similar  to  the  linear  interpolation  formul  (lo)  except 
that  the  factor  of  O.JOO  in  (ll)  is  13$  larger  than  the  factor  0.113  in  (.10),  so 
that  the  error  of  the  approximate  formula  (lo)  in  this  case  is  13$. 

For  the  case  that  the  outer  flow  is  a  rigid  rotation  witn  the  angular 
velocity  0,  i.e.,  v  =  rfi,  the  boundary  layer  Equations  (5)  have  been  solved 
by  Rogers  and  Lance  (i960)  for  several  values  of  —  . 
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According  to  tho  colution  ol'  Rogers  and  Lance  (i960),  the  radial  flow 
integral  la  given  by! 


00 

5  u  =  f  lid?.  =  .  ros  f(-) 

2  o  J  Vci)  'in' 


(15) 


where  the  function  f(^)  is  shown  in  Figure  2.  If  we  apply  our  approximation 
formula  (lQ)  to  the  case,  when  v  =  rO,  we  get: 


|-u6-=  '  ™  •  °’J^(1  -  |) 


Thus  the  function  f(— )  has  to  be  compared  with  the  approximate  expression  0.443 
(l  -  ^)j  which  ia, shown  by  the  dotted  line  in  Figure  2. 

The ’agreement  i a  still  good  enough  that  wc- can -consider  Equation  (lO)-as  a 
reasonable  upproximuliun  ulso  in  this  case.  Whether  or  not  Equation  (.10 )  ic 
approximately  valid  for  the  actual  velocity  profile’s  vo(r,t)  can  -be  checked 
after  obtaining  the  aolution  for  vn(r,t).  A  calculation  of  the  radial  flow 
integral  for  oome  of  the  obtained  velocity  profiles  vQ(r,i>)  has  been  done,  based 
on  the  simplified  momentum  integral  method  of  Mack  (1963)'.  The  uo  values  ob¬ 
tained  from  these  calculations  have  been  compared  with  the  approximate  valueB 
from  Equation  (10)  and  the  agreement  was  found  to  be  good  within  the  accuracy  of 
the  momentum  integral  method,  which  is  about  15$.  Thus,  we  can  conclude  this 
section  with  the  remark  that  the  error  of  the  approximation  for  uQ  (Equation  10) 

is  probably  not  larger  than  15$,  , 
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TT~  _ vi  h  -t  rtVi 

\J  EX.L1 (Jilt liV/  Ua  UJ-WU  AV*F 

written: 


a  ui 


-for  the  Reynolds  number,  .Equation  . (.10 )  can  be 


V  ‘  -  °'U3r  '^e  (rCD-  V 


(17) 


or,  with  the  notation 
k 

u_ 


k  (rcc  -  v  ) 


(18) 

(19) 


e.  Solution  for  v0(r,i.) 

After  substituting  (19)  into  Equation  (4)  the  equation  for  v  io: 

hv  v  dd  v  5v  Ir 

+  k(vo  -  +  ^ +  -it) 


l4 


?>v 

o 

8t 


(20) 


where  k  =  0. 443  jp  .  ^  . 

The  initial  and  boundary  conditions  are 

v  =  0  for  t  <  0 


v^  =  ato  for  r  =  a  and  t  >  0 
o  "  |; 

Let  us  restrict  our  analysis  at  this  instant  to  the  case  of  constant  cu 
(i.e.j  the  cylinder  is  started  at.t  =  0  to  spin  with  constant  angular  velocity 
o>).  In  many  cases  now,  we  can  neglect  the  friction  terms  at  the  right  hand  side 
of  Equation  ( SO-)  against  the  convection  terms.  To  see  this,  we  multiply 
Equation  (20)  with  — and,  introducing  the  dimensionless  variables  v* 


■  "  ani 

r*  =  —  Equation  (20)  becomes: 


any 


2 

"  4^  +  k(v*  -  r*)(|^  +  ^)  „  i_  (LZl  + 

5tbt  '  Sr*-  r*'  Re  '  ^2  + 


d  V*» 
"Sp* 


(21) 


with  the  boundary  conditions: 

v*  ■=  1  for  r* 


1  and  £bt  >  0 


It  can  be  seen  from  Equation  (2l)  that  the  solution  v*  la  a  function  of 
last,  r»  and  the  dimensionless  parameter  kRe  =  0.443  .  NRe  ,  i.e.,  v*  =  f(r*, 

tot,  kRe). 

If  kRe  ■  0.443.  r— ^Re  >  >  1  then  the  viscous  terms  at  the  right  hand  side 
,  .  n  V# 

of  (21)  become  small  except  for  small  times  when  the  gradient  is  large  near 

the  wall  r*  =  1,  11 

i  . .  « 

For  not  too  small  times  and  kRe  >  >1  wo  can  therefore  neglect  the  viscous 
terms  and  the  solution  of  the  inviscid  equation  is: 


r*e 


2tot 


v*  = 


1 

T'Jf 


2  tot 

e  -  1 


0 


for  r*  >  e 


for  r*  <e 


-tot 


•tot 


(22) 


A  plot  of  the  v*  profiles  (Equation  22)  for  different  times  is  shown  in 
Figure  3. 

It  is  remarkable  that  the  solution  (22)  satisfies  a(lso  the  complete 
Equation  (21)  with  viscous  terms  except  at  the  point  r*=  where  the  first 

derivative  is  discontinuous.  The  inviscid  equation  is  of  the  first  order  in  the 
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derivatives  and  hence  only  the  solution  function  has  to  be  continuous,  while  the 
complete  equation  has  second  order  derivatives  and  the  solution  must  have  con¬ 
tinuous  derivatives,  The  effect  of  the  viscosity,  therefore,  will  be  to  smooth 
-tot 


the  corner  at  r*  =  e 


With  the  solution  for  vq  (Equation  22)  the  other  flow 


components  uq. 


w  can  be  obtained  at  once, 
o 


First  of  all  we  have,  from  Equation 


(19): 


=  -  k(ro>  -  v  ). 


O  '  o' 

From  the  equation, of  continuity  (Equation  Id)  and  the  condition  that  the  flow" 
must  bc~  symmetric  to~the  plane~z  =  0,  it  follows:  ■  ■■  _ 

2ru 


w  -  -  — 


z~~  .0 

v  Sr  ' 


(23) 


Within  the  core  flow  we  can  distinguish  two  regions: 

,  -tot 
0  <  r/a  <  e 


-tot  .  „  /  .  , 

e  <  r/a  <  1 


Region  (l): 

Region  (2): 

According  to  Equation  (22)  the  particles  in  region  (l)  do  not  rotate  (i.e., 
v  o  0)  while  the  particles  in  region  (2)  rotate  with  the  velocity  vq  =  a or# 
given  by  Equation  (22).  Region  (l)  and  (2)  are  separated  by  the  cylindrical 
front  r/a  =  which  is  moving  toward  the  axiB  r  =  0. 

It  can  be  shown  that  the  particles  in  region  (l),  i.e,,  the  particles  ahead 

of  the  moving  front,  remain  ahead  of  it  until  they  hit  one  of  the  boundary  layers 
h 

at  the  cndwalls  z  =  +  ^  . 

To  see  this,  we  compute  the  trajectories  of  the  particles  in  region  (l). 

From  vQ  «  0  it  follows,  that  uq  -  -  kro>  and  wQ.  =-2kzu>  (Equation  23).  From: 


..  dr 
dt  "  uo 

dz 

dt  ”  wo 

we  obtain  by  integration: 


r  =  r  e 
o 


z  e 
o 


-  kreo 
2k  zw 

=-tot 

2tot 


(2)0 


/rO\2 

z  ( - ) 

err  ‘ 


(25) 


where  (rQ,zo)  is  the  particle  position  at  t  =  0. 

The  trajectories  (25)  course  entirely  in  region  (l),  so  that  the  solution 
(25)  is  compatible  with  the  supposition  v =  0. 


6 


\ 


Tims,  we  obtain  the  following  flow  picture: 

After  the  cylinder  is  started  at  t  =  0,  the  fluid  particles  move 
along  hyperbolas  given  by  (25)  until  they  hit  the  boundary  layer  at  one  of 

the  endwalls  z  =  +  —  . 

1.  —  ^ 

In  the  boundary  layer  the  flow  direction  changes  rapidly,  the  particles 
acquire  rotational  motion  and  are  thrown  radially  outwards  until  they  emerge 
from  the  boundary  layer  at  some  distance  behind  the  moving  front  r/a  = 
now  having  a  rotational  component  vQ  according  to  Equation  (22).  -  The  rest  of 
the  trajectory  courses  entirely  in  region  (2).  —  - 

Actually,  this  flow  picture  will  be  modified  slightly.  The  particles 
can  acquire  rotational  motion  already  in  region  (l)  by  the  action  of  the  viscous 
force  term  in  Equation  (21)  which  has  been  neglected  so  far.  This  will  also 
slightly  modify  ^he  trajectories- given  by  Equation  (25), 


f«  Equations  for  the  Angular  Momentum 

Of  particular  interest  is  the  loLul  angular  momentum  of  the  liquid  within 
the  cylinder,  which  In: 

ra  2 

I  »  fth.Sn  r  vQdr  (26) 

An  equation  for  the  angular  momentum  I  can  be  obtained  from  Equation  (20),  7 

To  thiB  end,  we  multiply  Equation  (20)  by  r  and  Integrate  from  r  ■  0  to 

Sv  v0  1  drv 

r  =  a.  If  we  oonsider  that  +  —  ■  -  ^  ,  this  glveB! 


;fe  l  r\ir  +  *  I  (rvo  - r  if  dr 


[v  f 


a  „  n  .  &rv 
2d  /I  _ o\. 

r  5  (F  "5F)dr 


(2*7) 


The  second  integral  consists  of  two  terms:  the  first  one  gives: 


ra 

J  rvo 


drv 

0 

dr 


dr»  = 


k  ,  f  >2 

2  (rV 


while  the  second  term  can  bs  integrated  by  parts: 


-  lao 


a  p  3rv 

r  '"y  dr 


-  kin  r  v 


f,u  2 

+  2kto  /  r  v  dr 

i  ° 
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Considering  that  vq  =  aw  for  r  =  a,  the  second  integral  of  Equation  (27)  gives: 


n  2  “"v 

k  /  (rvQ  -  r  w)  — dr 


k  4  2  „  pa  2  , 

-  r  a  a  +  2kw  /  r  v  dr 


The  integral  at  the  right  side  of  Equation.,  ( 27 )  becomes,  after  Integration  by 


parts; 


P  a  ^  t  drv  r  drv  -ia  a  *  &v„/r 

V  /  r2  I?  “ST)  dr  -  v  "  v(erv0]0  =  va3  (— gr) 

y.  -  o 


We  thus  get  the  equation: 


d  ra  2  ,  k 
dt J  r  vodr  "  2 


4  2  2  ’S  .5v0A> 

a  m  +  Eton  /  r  vQdr  =  vap  (— §—) 


According  to  (26)  the  value  of  the  angular  momentum  is  :  In  £ith 


'  2  . 
r  vQdr. 


In  the  final  state,  i.e.,  when  the  liquid  approaches  rigid  body  rotation 
v0  «  roi,  the  angular  momentum  becomes :  . 


T  =  2ith 

CO 


/a 


,  4 

r2wdr  «  2«h  .  p  . 


The  ratio  of  the  angular  momentum  I  to  its  final  value  *»>  then  is! 

I  I  r\dT 

7  °  “T7 - 

“  -a  m/4 

4 

After  dividing  Equation  (28)  by  it  can  be  written: 


§5  <f  >  *  t- 


i  J  E  4v 
J*-  J  aw 


itX  itd.) 

'  dr  ' 


If  we  neglect  again  the  viscous  term  at  the  right  side  of  Equation  (X>),  the 
equation  can  be  Integrated  at  once  to  give: 

l  -  (1  -  e’25^)  (3 

*  00 

Tills  result,  of  course,  could  have  djeen  obtained  also  from  the  solution 

i  v0 

Equation  (22)  for  the  velocity  profiles  v*  =  —  by  integration. 


irt 


In  order  to  improve  the  solution  (3l)  we  have  to  take  into  account  the 
.  Sv  /r 

friction  term  —  ( — «£-)  _  of  Equation  (JO).  This  can  he  done  within  a  fairly 

040  uT  T— 0 

good  approximation  by  assuming  the  approximate  shape  of  the  velocity  profiles 

dVjA  _  _ r _ 

vQ(r)  and  expressing  ( — ^r)r=a  l/l„- 

Since  the  contribution  of  the  viscous  tem  1b  small,  it  appears  reasonable 
to  assume  that  the  vq  profiles  have  essentially  the  same  shape  as  the  profiles, 
which  we  have  obtained  as  solution  for  the -inviscid- equation. - 

Thus,  we  may  assume  profiles  of  the  form: 


—  «  A  -  a/r  f 

A2  -  1 


aw 

vo 


^  =  0 

aw 


for  r/a  <  A 


-1 


(32) 


where  A  is  a  function  of  time.  For  the  case  that  the  viscous  term  1b  neglected 

V 

at  all,  the  solution  for  ^  was  given,  according  to  Equat(ion  (22),  by  the  profilea 
(32)  with  A  e  e1™^. 


It  might  be  mentioned,  that  the  vQ  profiles  (32)  satisfy  the  compatibility 

condition  p 

3"v.  ?4v  /r 

• .  (55) 


(' 


dr' 


i  ■*■*£> 


r=a 


0.. 


which  is  obtained  by  evaluating  Equation  (20)  at  r=a,  Using  the  vQ  profiles 
given  by  (32)  we  have: 


while 


AcA  _  £ks 

{— srW  n  “ 


1/1=1  -  A' 


A2  -  1 


Thus,  from  (31*)  and  (3?)  we  obtain: 

dv/r 
/  o  \ 

'  St  'r=a 


2u> 

a 


(34) 

(35) 

(36) 


Inserting  (36)  into  (30)  we  nave: 


ft  +  2to 

CO 


—  1  - 
I  , 

8v 

1 

k "  J 

~  2 
a 

m: '  \l 
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or,  after  dividing  by  to, 


il  -  8  r  1  / 

‘•k  r  *  1  Re  TJT  '  1 

u  CO  -  '  CO 


It  should  be  mentioned  that  the  approximation  for  (— jj? — )  ,  which  is 

given  by  Equation  (36),  is,, not  very  sensitive  to  the  special  choice  of  velocity 
profiles.  Although  the  very  early  v  profiles  are  somewhat  different  from  those 
fci-ven  in-  Equation-(32),  the  expression  (3.6)  is  still  a  good  approximation. 

It  can  be  obtained  from  Equation- (37)  that  for  sufficiently  small  times  the 
friction  term  is  dominant  however  large  kRe  is. 

The  equation  for  vQ  then  is: 


which  is  obtained  from  Equation  (£0)  by  neglecting  the  convection  .terms  and 


dv  /r 


compared  with  v 


The  solution  of  (38)  is: 


vQ(r,t)  -  ah) 


V-hf 


-y '  dy 


With  the  velocity  profiles  given  by  Equation  (39)  one  would  have: 


r*vn 


8  a)  l 

n  '  a  "  T/f 


2,54  ^  -iiz 


while  Equation  (36)  gives  approximately 


(¥0  - 


n  03'-  1 

2ai7i: 


Bince  for  small  times  we  can  neglect  1  against  jyj-  •  Thus,  we  see  that 
Equation  (36)  is  a  reasonably  good  approximation  even  for  the  very  early  velocity 
profiles. 

The  differential  Equation  (37)  now  has  to  be  solved  for  the  initial 
condition  i/l  =  0  at  t  =  0,  The  solution  for  i/l  can  be  given  implicitly: 
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Pkmt 


-1 


(to) 


1  + 


HT 

kRe 


log  (1  -  I/Ij  +  jgj-  log  (1  + 


kRe 


For  kRe  -> »  the  solution  (to)  reduces  to:  1  -  l/lM  =  e  i.e.,  the 

inviscid  solution  given  by  (Jl). 

For  very  small  times,  or  more  precisely  for  I <  <  1  and  i/l^  <  <  1, 

Equation  (to)  gives:  l/lw  K1  HJfte  •  ‘/tut;  i.e.,  the  angular  momentum  increases  as 
the  square  root  of  t\., _ •  • _  _  _ _ 


The  validity  of  the  preceding  results -is-reatricted  to  the-case  where  the 
Angular  velocity  of  the  cylinder,  oi,  remains  constant  after  the  cylinder  Is 
Itarted, 

If  oi  i3  not  constant,  it  is  advantageous  .to  introduce  the  dimensionless 
quantity 


!*  =  ± 


a  oi 


SJ 

vor,  dr, 


\\  - 


TO 


where  aiQ  is  a  constant  reference  angular  velocity  and  Iq  is  a  reference  angular 
momentum : 


2jt  p 


h  f  %r3<ir' 


the  latter  corresponding  to  a  rigid  rotation  with  angular  velocity  u>  . 


r 


Dividing  Equation  (28)  by  — r —  and  substi 
to  Equation  (36)  we  have; 


“V - 


according 


r«a 


dl* 

a2a> 


t'2k  ( 


01 


o  '01, 


A  r, 


I*  .  2L 

8 

01 

L  ai  , 

rrr  ““  ■  ± 

“o. 

"  Reo 

“o 

I*  CD 

0 

(to) 


where:  Rc\_ 


and  k0  -  O.MO  ^  .  (Re„) 


rl/2 


Equation  (to)  has  been  used  to  calculate  the  spin  decay  of  a  liquid-filled 
shell.  If  a  shell  containing  a  liquid-filled  cylindrical  cavity  is  started  im¬ 
pulsively  to  spin  about  its  axis  of  rotation,  the  liquid  continuously  absorbs 
angular  momentum,  thus  reducing  the  spin  until  the  liquid  finally  attains  rigid 

I 

body  rotation, 
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Let  A  be  the  axial  moment  of  inertia  of  the  empty  shell,  o>  the  instantane¬ 
ous  angular  velocity  and  I  the  angular  momentum  of  the  liquid,  then,  due  to  the 
conservation  of  angular  momentum, 

I  +  txA:  =  constant  =  cdqA,  ( 43  ) 


where  idq  is  the  initial  angular  velocity  when  the  liquid  is  at  rest. 

According  to  (1)3)  we  can  express  1  by  cut 
I  =  Acuq(1  -  o/o>0), 
or  using  the  definition  (4l): 

Acs  _ 

I*  =  (1  -  cu/cu0).  (44) 

!  • 

After  inserting  (44)  intv)  (42)  we  get  a  differential  equation  for"l*  (or  <u/cuQ) 
which  may  bo  solved,  by  numerical ..integration. 

g,  .Turbuletft. Boundary  Layer  Flow  ■  -  ..  ______  •  . 

In  Sections  (c)  and, (d)  we  have  assumed  that  the  boundary  layer  flow  at  thg 

cylinder  ends  is  laminar.  This  assumption  is  valid  for  Reynolds  numbers  Re  =  S-iS 

ft  \>  v 

up  to  about  3  ‘  103  (see  e.g.  Schliehting  1930)-  For  Reynolds  numbers  greater 

than  3  ■  10^  the  boundary  layer  flow  will  be  turbulent  and  the  radial  ma3s  flow 

and  hence  uQ  will  be  different  from  the  laminar  case. 

According  to  the  solution  of  von  Karm,an  (1921)  for  the  turbulent  boundary 
layer  flow  on  a  rotating  disk,  the  radial  flow  integral  is: 


This  formula  is  analogous  to  Equation  (_6) >  which  corresponds  to  the  laminar  case. 

;j  g  | 

Using  again  the  notation  Re  -  for  the  Reynolds  number,  Equation  (45) 
can  be  wri'', cent 


By  arguments  similar  to  those  used  in  Section  (d)  for  the  laminar  boundary 
layer,  it  seems  appropriate  to  generalize  Equation  (46)  to: 


.  /f-  (rcu  -  v  )  ' 

0.035  a  .  (Re)-1^ - tvr- 

(acn  )5/5 


( **7 ) 
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According  to  Equation  (6)  we  then  have  for  the  radial  component  of  the  core  flow: 

bb 


n  ,  /l=  (ra>  -  v  ) 

u  (r)  =  -  0.035  •  (tt)  .  (Re)'1/5  ° 


(«n) 


3/5 


Equation  (48)  is  analogous  to  Equation  (17)*  With  u  obtained  from  Equation  (48) 
the  equation  for  the  core  flow  (Equation  4)  becomes: 


-dv 

o 

IT 


!.  - q/c-8v--  v 

J75  (r“"  vo}  '(“5F  +  r] 


where : 


(ad.) 


0.035 '(^)  (Re)-1/5 


Svo/r-, 


(49) 


(50) 


Analogous  to  the  procedure  described  in  Section  (f),  we  find  the  equation  for 
the  angular  momentum  by  multiplying  Equation  (49)  with  r^  .  — r— -  and 


r? 

a“o 


integrating  over  r. 


If  we  further  introduce  again  the  dimensionless  quantity  I*  according  to 
Equation  (4l)  we  finally  have: 


dl* 

dib0t 


*t  4 

(aiu)5/5  a^to  S 


£  r(*D  -  v//5(*g)  dr'.  (51) 


The  integral  in  Equation  (51)  cannot  be  evaluated  as  it  could  for  the  laminar 
case,  without  knowing  how  vQ  depends  on  r.  Only  for  small  times  when  the  v0 
profiles  are  restricted  to  a  nanwrzone  near  the  wall  r=a,  we  can  approximate 
r  in  the  Integral  by  a  and  find: 


/ 


a  0  /_  5a  v 

*<•»  -  *„>  3/5  -y-  ^  - 


j|  «2  («  -  «,)l5/5 


?J  .3-)13/5 


Inserting  the  last  result  into  (51)  and  using  for  the  viscous  term  at  the  right 
side  the  same  approximation  as  in  Equation  (42),  we  obtain:  1 


'j£L  .  k  (^-)9/5  =»  r-V  (-) 

cG£t  T 5  Kto  'o>0'  (Re)0 


1  u>  , 

I*  «b  ’ 


(52) 


where  again  (Re)Q  =  and  kfco  =  0.035  (|p)  (Re0)-1/5 

In  order  to  evaluate  the  Integral  in  Equation  (51)  for  later  times,  when  the 
assumption  r  «  a  is  no  longer  valid,  we  have  to  make  assumptions  about  the  shape 
of  the  velocity  profiles  v0(r).  We  may  assume,  that  the  velocity  profiles  are 


given  again  roughly  by  Equation  (32).  These  profiles,  at  least,  satisfy 
the  compatibility  condition  (33).  i.e.,  they  are  correct  near  the  wall. 

Using  the  vQ  profiles  given  by  (32)  the  integral  in  Equation  (51)  can  be 
evaluated  and  expressed  by  I*.  After  some  lengthy  calculation,  one  finally 
obtains  from  Equation  (51): 


ai*  .  k  f£Lo9/5 

druot 


^8/5 


O 


9_  to_  1  co 
so  “o  -1*  “o 


For  small  times,  i.e.,  as  long  as  I*  —  <  <  1,  the  integral  approximates  to: 

k  u>  13/5 

(I*  — )  , 

13  CD  _  '  _ 

‘■3  %  8/5  4 

while  (1  -  I*‘ jp)  «  1,  and  Equation  (53)  reduces  to  Equation  (52). 

Equation  (53)  has  been  used  to  calculate  the  spin  decay  of  liquid-filled 

shell  for  Re  >  3.1o5. 
o 


III.  COMPARISON  WITH  EXPERIMENTS 

In  order  to  test  the  analysis  given  in  the  preceding  section,  some  of  the 
theoretical  predictions  have  been  compared  with  existing  experimental  data. 

A  detailed  description  of  the  experimental  arrangements  is  given  by 
Karpov  (1962), 

A  quantity  which  has  been  measured  directly  is  the  axial  spin  decay  of 
liquld-f tiled  shells,  fired  from  a  gun.  After  the  shell  leaves  the  gun,  the 
angular  velocity  decreases  continuously.  The  decrease  of  angular  velocity  Is 
caused  by  absorption  of  angular  momentum  in  the  liquid  and  by  the  torque  due  to 
air  friction.  The  contribution  of  the  air  friction,  which  is  usually  small,  can, 
be  determined  separately  by  observation  of  the  spin  decrease  of  the  empty  shell. 
Tine  difference,  -which  is  due  to  absorption  of  angular  momentum  has  been  plotted 
for  two  typical  cases  in  Figure  1  and  Figure  5.  For  comparison,  the  theoretical 
curves  are  plotted  in  the  same  diagram  and  also  the  curves  obtained  from  the 
theory  without  secondary  flow.  The  fineness  ratio  of  the  cylindrical  cavity  was 
in  both  eases  =  2.68.  For  fineness  ratios  smaller  than  2.68  the  effect  of  the 
r.econdury  flow  would  be  even  more  pronounced. 
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The  Reynolds  number  for  the  ease  plotted  in  Figure  U  was  Re  =  I.76  x  10'3  so 
that  a  laminar  boundary  layer  flow  could  be  assumed,  while  for  the.  case  of 
Figure  5  the  Reynold's  number  of  Re  =  6.1  x  lCr  was  above  critical  and  therefore 
the ;i  formula  for  turbulent  boundary  layer  flow  was  applied. 

In  addition  to  observations  of  spin  decay,  .experiments  have  been  done  to 
observe  the  secondary  flow  itself.  To  this  end  an  impulsive  spin  generator  was  .. 
designed  which  is  described  by  Stoller  (i960).  ;The  spin  generator  consisted  of 
a  liquid-filled  cylinder  with  transparent  walls  .which  could  be  started  impulsively 
to  spin  about  its  axis.  A  suspension  of  small  particles  was  dissolved  in  the 
liquid,  the  specific  gravity  of  the  particles  being  the  same  as  that  of  the 
liquid.  The  trajectories  of  the  particles  could  be  observed  with  the  aid  of  a 
motion  camera.  For  this  type  of  observation,  of  course,  only  transparent  liquids 
could  be  used. 

Figure  6  shows  some  of  the  observed  parLicle  trajectories,  or  rather  the 
first  part  of  it,  and  the  particle  positions  at  constant  time  intervals.  Close 
to  the  observed  trajectories,  theoretical  trajectories  and  particle  positions- - 
according  to  Equation  (25)  are  plotted.  The  agreement  is  reasonably  good. 

Some  deviations  can  be  explained  by  the  fact  that  the  Reynolds  number  of  the 
experimental  flow  was  rather  low  (Re  =  1,83  x  10^)  while  the  .theoretical 
prediction  is  based  on  the  assumption  of  high  Reynolds  numbers,  where  the  viscous 
forces  in  the  core  flow  can  be  neglected. 
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